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ABSTRACT 

The axial anomaly equation in 1+1 dimensional QED is obtained on the lattice for fermions obeying 
the Ginsparg- Wilson relation. We make use of the properties of the Lattice Dirac sea to investigate the 
connection between the anomaly and the Ginsparg- Wilson operator in the Hamiltonian picture. The correct 
anomaly is reproduced for gauge fields whose characteristic time is much larger than the lattice spacing, 
which is the regime where the adiabatic approximation applies. A non-zero Wilson r parameter is necessary 
to get the correct anomaly. The anomaly is shown to be independent of r for r > 0.5. The generalization 
to 3+1 dimensions is also discussed. 
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The lattice regularization is one of the few non-perturbative methods available for defining quantum 
field theories. Lattice gauge theories have revealed many interesting features of gauge theories that are not 
easily visible in the usual perturbative approach. Nevertheless, the lattice regulator has proved problematic 
if fermions have to be incorporated into the theory. A naive discretization of the fermionic theory suffers 
from the replication of fermion modes due to the " doublers" . The doublers are degenerate in energy with 
the originally introduced fermions and though they have lattice momenta of the order of the cut-off (1/a in 
lattice theories, a is the lattice spacing), they mimic ordinary low energy fermions. These doubler modes 
cannot be ignored as they participate in physical processes, for instance they can be pair created, and can 
affect the value of physical quantities- such as the free energy. The first method to handle these doublers 
was given in [Q and it uses an additional term in the action -the Wilson term- to lift the degeneracy 
of the fermions, thereby decoupling the doublers in the continuum limit. However, this method has the 
disadvantage of explicitly breaking chiral symmetry, and hinders the study of dynamical questions related 
to chiral symmetry breaking. A cure for the doubling problem that explicitly breaks chiral symmetry also 
makes the lattice regularization of chiral gauge theories, such as the standard model, much more difficult. 
The Nielsen-Ninomiya no-go theorem [|j decrees that any chirally symmetric lattice Hamiltonian satisfying 
general properties like locality and hermiticity must result in a replication of fermion species. This theorem 
seems to suggest the impossibility of defining undoubled fermions on the lattice without breaking chiral 
symmetry. Recently, however, alternative methods for tackling these problems have emerged. One of them 
1^ uses the so called Ginsparg- Wilson relation [Q for Dirac fermions. The Ginsparg- Wilson (G-W) operator 
is obtained by the application of block-spin transformations to a chirally invariant Dirac operator, and 
which therefore suffers from fermion doubling, using a chirally non-invariant blocking kernel. Although 
the G-W operator is not chirally invariant, it contains the information of chiral symmetry because it has 
been obtained after blocking a chirally invariant lattice action. Its construction, by a renormalization 
group transformation of a chirally invariant action, is bound to leave the low energy properties related 
to the chiral symmetry unchanged. This approach of formulating lattice fermions has led to many recent 
developments such as, lattice formulations of chiral symmetry, the search for chiral gauge theories on 
the lattice, methods of defining a lattice topological charge, and formulation of lattice index theorems etc. 
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The G-W operator has to satisy the following relation 

^75 +75^ = 0^75^5 ■ (1) 

(There are different versions of the G-W relation depending on the precise form of the blocking kernel used. 
The above form is one of the simpler ones and is sufficient for the ensuing discussion. Here a is the lattice 
spacing. ) The G-W operator clearly does not satisfy chiral symmetry (because {D, 75} ^ 0). Even though 
the G-W operator seems to share the properties of a chirally noninvariant mass term, it is a milder way to 
break the chiral symmetry on the lattice. This is because it is obtained by blocking a chirally symmetric 
action, although using a chirally non-invariant kernel. The low energy properties of the G-W operator on 
the lattice are the same as those of the chirally symmetric action. 



Another approach to problems of chirality on the lattice is the overlap approach introduced in |15|. It 
captures many essential elements of domain wall fermions as well as the one in [ p^ which requires 
an infinite number of auxiliary fields. The chiral determinant is expressed as an overlap of the ground 
states of two many body Hamiltonians and a construction of chiral gauge theories involves regularizing 
the overlap . Though the original Ginsparg- Wilson approach and the overlap approach appear to have 
nothing in common, the overlap operator (which appears in the Hamiltonian) has been shown to satisfy 
the Ginsparg- Wilson relation ||l9|| . The overlap is also a real time approach since it involves the quantum 
mechanical scalar product of the ground states of two Hamiltonians. This approach has led to many further 
studies of chirality on the lattice |Q . 

Another useful way of looking at the fermion doubling problem on the lattice is to look at the chiral 
anomaly structure of the lattice theory and to see what it yields in the continuum limit. A symmetry is 
said to be anomalous if it is no longer present in the quantum theory although it is present in the classical 
theory. The anomaly manifests itself by a non-conservation of a classically conserved charge. Anomalies 
are an inescapable part of some quantum field theories and have many important physical consequences. 
Their origin is related to the problem of regularizing amplitudes in quantum field theories while maintaining 
their invariances. In the path integral formulation of quantum field theory they arise because of the non- 
invariancc of the measure of the path integral As is well known 1+1 dimensional QED has a chiral 
anomaly when massless fermions are present. The anomaly arises because it is not possible to find a 



regularization of the gauge theory which maintains both the gauge invariance and chiral symmetry. A 
simple way of demonstrating the anomaly in 1+1 dimensions is by using a gauge invariant point-split 
definition of the axial vector current 0] which can be seen not to be conserved. There are no anomalies 
in the naive latticisation because it is a gauge invariant regulator which also maintains chiral symmetry, 
but it is impossible to put only fields of one chirality on the lattice. This is consistent with the fact that 
no regularisation exists which simultaneously preserves gauge invariance and chiral symmetry for arbitrary 
matter content. The continuum limit of the lattice theory, on the other hand, must be able to reproduce 
the correct anomaly structure of a given theory. The naive fermionic lattice action coupled to gauge fields 
gives an anomaly free theory in the continuum limit because the anomalies are cancelled between the naive 
and doubled modes [|ll| . It was shown in ||ll[ that the Wilson term reproduces the correct anomaly on the 
lattice provided the symmetry breaking parameter r 7^ 0, the anomaly in the continuum limit being given 
by the co-efficient of the Wilson term. This is quite a surprising result because the Wilson term explicitly 
breaks the chiral symmetry but yet reproduces the anomaly which is essentially a quantum mechanical 
breakdown of the classical chiral symmetry. 

As stressed by Nielsen and Ninomiya , and Peskin, the Hamiltonian formulation provides a much clearer 
physical picture of the anomaly in terms of the energy level shifting of the filled Dirac sea ( for a very 
clear exposition see also p^). In this picture the anomaly arises because pairs of net chirality are pumped 
out of the infinitely filled Dirac sea. If one tries to transcribe this picture on the lattice, as was done 
by Ambjorn et al |l^, one finds that the lattice Dirac sea is always finite and the anomaly always gets 
cancelled by the doubler modes in the absence of the Wilson term [^l]jlj]. The physical picture of the 
anomaly presented in |l2j can be applied on the lattice with a Wilson mass term. The role of the Wilson 
mass term is to suppress the contributions to the chiral charge coming from the doubler modes resulting 
in a non-zero anomaly on the lattice plf . 

Our aim is to carry out a similar analysis for fermions satisfying the G-W relation and to see how a 
non-zero anomaly comes about on the lattice. This should complement the derivation of the anomaly from 
the Ginsparg- Wilson action in the Euclideanised formalism where the anomaly is showed to arise out of 
the measure B. We mention here that the axial anomaly in 1+1 dimensional QED is also reproduced in 
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the overlap formulation | [l5| . Our derivation, apart from being quite different from the methods employed 
in , also highlights the role played by the Wilson term in giving the correct anomaly. The discussion 
will be in the Hamiltonian framework and we will derive the anomaly equation for the abclian theory in 
the 1+1 dimensions. We will then comment on the extension of this picture to 3+1 dimensions. 
The Lagrangian density for a Dirac fermion in Minkowski space is given by 

L{ip, -0) = -0(^7^.5^ - m)ip ; (2) 

the gamma matrices satisfy 7q — 70, 7I = ~7j:7 and obey the relation {7^,7,/} — ~2g^y. The metric g^y 
is diag{l, —1, —1, —1). ip{x) denotes the relativistic adjoint ip''{x)jo. We shall use the Weyl representation 
for the gamma matrices. The Hamiltonian density is given by 

H = iP{x){-nA + m)ip{x) (3) 

Before we discuss the GW fermions in the Hamiltonian picture, it is instructive to briefly review how 
the unwanted doublers are handled in the Euclidean formalism with the help of the Wilson mass term. In 
Euclidean space the Lagrangian density becomes 

= ij{-ij^df, - m)tp] (4) 

the Euclidean gamma matrices satisfy {7^,71/} = —25^1, , and 7^ = — 7/j. The naive lattice discretization 
of the Dirac action for a massless fermion in Euclidean space is given by 

S = ^-^i^{x)-it{^{x + i)-'ip{x-i)) . (5) 

X i 

In momentum space (in d dimensions) this becomes 

I d''fc^(fc)(^(5i^h,)^(A) . (6) 
Jbz a 

BZ denotes the range of integration to be the d dimensional Brillouin zone of the lattice. As is well known, 
the above discretization suffers from the presence of additional fermions at the corners of the Brillouin zone 
(2** in d Euclidean dimensions) leading to 2'^ fermions in the a limit. A method for eliminating the 
unwanted fermions is to give them very high masses in the continuum limit. The oldest way of achieving 
this is by adding a Wilson term to the massless action 



- — ^ ij{x){'ijj{x + i)+ iPix 2'ijj{x)) 



(7) 



(The Wilson term mimicks a mass term although in a more subtle way; the mass terms are "momentum 
dependent") This leads to the modified propagator (in momentum space) 



In the above expression we have also introduced a bare mass term m. The momentum dependent mass 
terms ensure that the modes at the corners of the Brillouin zone have masses of the order of 1/a and 
decouple from the low energy effects (in the limit a — s- 0). The price paid for eliminating these doublers is 
the lack of chiral symmetry in the fermion action (the Wilson term explicitly breaks chiral symmetry ). 

In order to define the real time evolution of Ginsparg- Wilson fcrmions in the Hamiltonian formulation, 
we have to first construct a Hamiltonian operator starting from the Euclidean functional integral. In the 
transfer matrix formalism this is done by choosing a particular axis as the time direction (with a lattice 
spacing r) and then taking the so called r continuum limit. The r continuum limit (r 0) is taken on an 
anisotropic lattice with different spacings in the space and the time directions. In the Euclidean formulation, 
the Ginsparg- Wilson operator can be formally understood to have arisen out of a chirally non-invariant 
blocking transformation. One could have used a r-continuum procedure to obtain the Hamiltonian version 
of the GW prescription. However, since we are eventually interested only in the Hamiltonian, it is sufficient 
to do the block spinning only in the spatial directions which already yields an (anisotropic) lattice with a 
blocked action. This blocked action has the same partition function as the original action, though spatial 
correlation lengths are halved on this lattice. We can then proceed to construct a Hamiltonian operator 
via the transfer matrix in the usual way for fcrmions This way we are able to study the real time 
evolution of Ginsparg- Wilson fcrmions. As is well known, the passage from the transfer-matrix formalism 
to the Hamiltonian formalism can be carried out in more than one way, depending on the simplicity or 
complexity one desires. Our Hamiltonian is invariant under the symmetry discussed by Luscher in Q and 
is a valid starting point for carrying out calculations of wave functions, ground states, excited states, etc. 
The chirally invariant action on the lattice is 
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d 




(8) 




(9) 



m,n 
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where h satisfies {h, 75} = 0. We can block only the spatial degrees of freedom by using the kernel (defined 
in|]) 

K{ip,(p) = exp-^^{ij{m) - x{m)){ij{m) ~ xim)) ■ (10) 

mn 

The only difference is that the blocked field x(to) is defined as 

X(m,t) =77^(/.(n,i) , (11) 

n 

where the above summation is only over fields defined over a hypercube at the same instant. This is the 
condition which ensures that only the spatial degrees of freedom are blocked. This leads to the following 
effective action for the fermions on the blocked lattice 



exp{-A[iPi;]) ^ J d4>d(j)eyip[-A{(l3,4>)]K{i;,(j)) . (12) 
The action on the blocked lattice is given by 

=^V;(m)/i:„„V(n) , (13) 

m,n 

where the fields V' are defined on the blocked lattice which has twice the spacing (in the spatial direction) 
of the original lattice. The propagator of the blocked lattice satisfies the following relations 

{ht,l5} = 

After the usual passage to the Hamiltonian via a transfer matrix formalism the Hamiltonian one obtains 
is simply 

H = J2^\^hohs{x,y)^{y) . (14) 

Henceforth, the operator hg will be called D and it satisfies the G-W relation. Any hg satisfying the 
Ginsparg- Wilson relation can be used for defining our Hamiltonian. An explicit choice for D satisfying the 
Ginsparg- Wilson relation is the overlap operator given by Neuberger's construction Il5| 



A is defined in terms of the Wilson operator Dyj as 



A=l-aD, 



W 



(16) 



It can be easily checked that the operator D satisfies the G-W relation. It is worth adding here that any 
A satisfying the following properties: = 75A75 and A^A commutes with A will give a D satisfying the 
GW relation. This may be useful in a more general context. 

Though any D satisfying the Ginsparg- Wilson relation can be used to construct our Hamiltonian, we 
have used the above explicit form proposed by Neuberger et al for our calculations. It should be stressed 
that apart from this choice, the considerations of this paper are independent of the overlap formalism. 
The above operator relations can be translated into momentum space and it is in momentum space that 
we will make most of our manipulations. It should be emphasised that in general, with external fields, 
momentum space description is not very economical.But in the 1 + 1 dimensional abelian case with only 
an electric field, and in the 3 + 1 dimensional case with uniform electric and magnetic fields, momentum 
space description is still useful. 

The Hamiltonian of the lattice field theory is 



x,y 

•ijj and V can be interpreted as field operators in the usual sense. The wave equation for the fermion fields 
is 



Using the properties of D and 70 it is easy to show that the Hamiltonian is hermitian, and therefore 
the evolution is unitary. We are basically interested in how the anomaly arises in this model. To get the 
anomaly we must of course couple the fermions to an external gauge field and then look for non-conservation 
of the chiral charge. It is well known that the anomaly can be extracted by treating the gauge fields as a 
classical variable and quantizing only the fermions. In order to be able to study the problem of fermions in 
an external field we will have to make some approximations which will be described shortly. The analysis 
will be presented in 1+1 dimensions. 




(17) 




(18) 



y 
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In 1+1 dimensions the only effect of the external gauge field is to shift the momentum variable k to 
ga{t) where a{t) is the time dependent component of the vector potential (yli(x,t) — a{t) andAo — 0). 
This is true only in the Hamiltonian picture, even in the 1 + 1 case, because in the Euclidean case, 
D'' D = (fc — ea(t))^ + 75^, where £ is the electric field. We define the chiral charge operator on the lattice 
in the usual manner as 

Q„ = a ^V^(x,t)75^(x,t) . (19) 

X 

All the operators are defined in the Heisenberg representation. The time dependence of the chiral charge 
is given by 

Q„ = ^+*[ff,Qn] . (20) 

ot 

Since only the second term on the righthand side contributes to Qn we have 

Qn^t[H,Qn] . (21) 

Evaluating the commutator this becomes 

Q„=iY,i^ix,t){D,l5mx,t) . (22) 
Using the Ginsparg- Wilson relation this becomes 

Q„ = za ^ {^{x, t)Dj5D^{y, t) . (23) 

x,y 

The spatial indices of D and all Dirac indices have been suppressed for ease in reading. Since Dj^D is 
non-zero we see that the chiral charge is not in general conserved, as expected. It remains to be seen if 
this non-conservation of the chiral charge reproduces the correct anomaly. For this the Dirac equation in 
Eq. ^ has to be solved in the presence of an external potential and the solutions have to be examined. 
We consider a time dependent potential which rises from zero to a constant value Ar in a time r. r is a 
time scale in the problem and two cases can be easily analyzed, the sudden limit r — s- and the adiabatic 
limit T — !■ oo, as was also done in p^ . First one writes the free field ^p{x,t) as a superposition of positive 
and negative energy spinors 

ip{x,t)^[ ^[b{k)u{k)exp{~ikx) + d\k)v{k)exp{ikx)] . (24) 
Jbz 27r 



kx is short for koxo — kixi. E = fco- As usual u{k) and v{k) represent positive and negative energy 
spinors. The operator b{k) destroys an electron of momentum k\ and the operator d^k) creates a positron 
of momentum —ki. Putting these spinors in the equation for the axial charge we have 

{\Qn\)= ! dk^v{k,t)D{k,t)^^D{k,t)v{k,t){d{k)d)(k)) . (25) 

Jbz 

The angular brackets denote expectation values in the vacuum which is the state with zero electrons and 
positrons. Using the definition of D and the properties of the gamma matrices it is easy to show that 

D{k)-i^D{k) = ^^D\k)D{k) . (26) 

Now , D'^ {k)D(k) is a c-number and acts trivially on the Dirac spinors. {D'^D is not a c-number in 
general. For the specific class where D is of the form Ai^yi + B with Ai, B commuting, this is true. Already 
in the 3 + 1 dimensional case this is no longer true even for a uniform magnetic field.) In order to evaluate 
v{k, t)j5v(k, t) we have to determine the evolution of the negative energy spinor in the external field a{t). 
Before we calculate this quantity it is instructive to calculate the same without any field. In the absence 
of an external field the positive and negative energy states evolve as 

u{k,t) = exp{—iE{k)t)u{k) 
v{k,t) = exp{iE{k)t)v{k) . 

The spinors u{k) and v{k) satisfy the time independent Schroedinger equations with positive and negative 
energies. 

joD{k)u{k) = E{k)u{k) 
-foD{-k)v{k) = -E{k)v{k) . 

The Weyl representation for the 2 dimensional 7 matrices is 

71 = h5 = io-3 70 = fi 75 = (72 . (27) 

(Ti and 0-3 are the Pauli matrices. The eigen values of the spinors are given by E{k)'^ = D\k)D{k). Using 
the definition of D{k) we get D^{k)D{k) to be a c-number. D{k) can be written as 
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If we now write D{k) as 

D{k)=g{k,t)+jif{k,t) , (29) 
then E'^{k) = D\k)D{k) = P{k,t) + g'^{k,t). f and g are functions given by 

,(.) = i(l-ii^^) 

a^^{A^k)y 

The function / is an odd function of k whereas the function g is an even function of k. The time independent 
spinors can be normahzed to satisiy 

u{k)'y5u{k) = 
v{k)j5v{k) = 
u{k)j5v{k) = 1 . 

This means that (Qn) = in the absence of an external field. This is as expected, there is no anomaly in 
zero external field even though the Hamiltonian is not chirally invariant. 

In the presence of an external field the only change in the operator D is a replacement of by fc — 17 a(<:). 
Since the structure of D is not afi'ected by an external field the Ginsparg- Wilson relation is still satisfied. 
Although the evolution of Dirac spinors in an arbitrary external field can only be analyzed numerically, 
two limiting cases admit a simpler analysis. These are the adiabatic limit and the sudden limit, and we 
shall examine these two cases separately. When an external field is turned on slowly (compared to the time 
scales in the system) we can use the adiabatic approximation. The sudden approximation is useful when 
the field is turned on faster than the fastest time scale in the system. The rate at which the field is turned 
on can be controlled by introducing a parameter r defined as follows 

A{t) = i < (30) 



A{t) = A t> T 
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(31) 



The precise form of A(t) for < i < r is not very important. The sudden Umit corresponds to t ^ and 
the adiabatic hmit corresponds to r ^ cx3. In the adiabatic approximation the form of the positive and 
negative energy spinors for times < i < r is given by 

u(k, t) = a{k, t)C{t)u°{k, t) + b{k, t)^{t)v°{k, t) (32) 

v{k,t)^c{k,i)C{t)u°{k,t) + d{k,t)£,{t)v°{k,t) . (33) 

The above equation is written for the individual fourier components of the positive and negative energy 
spinors in an external field, v?{k,t) and v'^{k,t) are the positive and negative energy spinors satisfied by 
the instanteneous Schrodinger equation at the instant t. We closely follow the notation of and we 
have also corrected some of the misprints which occur therein. S,{t)* and S_{t) are the phase factors for the 
positive and negative energy spinors. The phase factor S^it) is given by 

£,{t) = exp(i / E{k, t')dt' . (34) 
Jo 

a,b,c, and d are called Bogoulobov coefficients. The Bogoulobov co-efficients are chosen to satisfy the 
following boundary conditions 

a(0) = d(0) = 1 6(0) = c(0) = . (35) 

These boundary conditions ensure that we are looking at the evolution of the positive and the negative 
energy states before the external field is switched on. As mentioned before, u°(fc, t) and u°(fc, t) are positive 
and negative energy spinors having momentum k and —k respectively, and they satisfy the instanteneous 
Schrodinger equations given by 

-foD{k)u°{k, t) = E{k, t)u°{k, t) 
7o£'(-fc)w°(fc, t) = -E{k, t)v°{k, t) . 

Substituting the expressions in Eq. ^ and Eq. |3^ in the wave equation for the fcrmions and using the 
previously mentioned boundary conditions we get 

c{k,t)^ f a{k,t')d{k,t')S,'^(t')dt' (36) 
12 



t 

(i(fc,t) = l-/ a{k,t')c{k,t')£,*'^{t')dt' (37) 
Jo 

b{k,t) = -c*{k,t) (38) 
a{k,t) = d*{k,t) . (39) 

The quantity a{k,t) is defined by 

u°{k,t) = a{k,t)v°{k,t) , (40) 

and is 

After using the stated normalizations of the spinors, v{k,t)j^v{k,t) is given 

v{k, t)j5v{k, t) = cdi'^ - c.c . (42) 

The co-efScients c{k,t) and d{k,t) can be approximated by (for small values of a(k,t)) 

c(fc,t) = -*^|^^'(0 d{k,t) = l + Oia') . (43) 

In the adiabatic approximation the quantities inside the integrand on the right hand side of Eq. ^ and 
Eq. ^ are evaluated at t — 0. Substituting for the values of /(fc, t) and g{k, t) and using the relation 

Dik)jr,Dik) = j5D\k)Dik) = ^5E\k) (44) 

the r.h.s of Eq. ^ becomes 

/ ^dk C{k) (a a ga{t)) (45) 
Jsz 27r 

where C(fc), a complicated expression, is given in the appendix along with the expressions for /, g, /, g. The 
function C{k) can be plotted and the integral of C(fc) over the Brillouin zone can be estimated numerically. 
The function C{k) depends on r,k,a,m. We plot C{k) as a function of k in Fig. to Fig. . We first plot 
it for zero mass and then for a non-zero value (m = 5). The eigenvalues of the spinors (in a zero external 
field) are also plotted in Fig. to Fig. . The first thing we observe is that when r = the eigenvalue 
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spectrum docs not distinguish between the modes at fc = and k = tt/o, and the integral of C{k) over the 
Brillouin zone is zero. For r =/= the modes at fc = and k ^ ir/a have different energies and an asymmetry 
develops in the function C{k). For r = 1, the integral of C(fc) over the Brillouin Zone has the value —2. 
Substituting this in Eq. we get the anomaly equation 



in the continuum limit. 

We have studied the function C{k) for different values of r and we find that the anomaly is independent 
of r for large r but vanishes for a smaller and, in particular, a zero value of r. The value r = 0.5 seems 
to separate the region with and without the anomaly. A non-zero Wilson r parameter is necessary to get 
the anomaly in the continuum limit. This means that a Neuberger like operator for D where the naive 
Dirac operator is used in place of Z?„, in Eq. |l6| will not reproduce the correct anomaly in the continuum 
limit inspite of satisfying the Ginsparg- Wilson relation (the operator D with D^; replaced by Dnaive also 
satisfies the Ginsparg- Wilson relation). The case m 7^ can also be analyzed along the same lines and it 
turns out that the integral of C(fc) is very small, consistent with zero. It appears that in this case we have 
an exact cancellation of the bare mass term with the anomaly term to give a zero rate of change of chiral 
charge. Nevertheless, the anomaly term is still present and so is the mass term, but the two appear with 
opposite signs. The adiabatic approximation is justified when the switching time r is much greater than 
the characteristic time periods of the system. In our example 2tt / E{k) is the characteristic time period of 
the system and the adiabatic approximation is justified when t >> a. To summarize, a zero bare mass 
term with a non-zero r parameter gives an anomaly independent of r for r > 0.5. When a bare mass term 
is included we have a cancellation of the anomaly term with the mass term though both terms are still 
present. 

The rate of change of the chiral charge can also be calculated in the sudden approximation in which an 
external field is turned on infinitely fast. This approximation corresponds to the limit r ^ 0. In this limit 
the spinors vP{k, t), t) arc unchanged immediately after the field is turned on and only their evolution 
is governed by the new Hamiltonian (with a constant field) . The normalization of the spinors in a constant 
external field (with value A^) can be made just as in the zero field case and there is no anomaly in this 




(46) 
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limit. This limit corresponds to the case t << a. 

The time r is a characteristic time associated with the gauge fields on the lattice and our calculation 
clearly shows that in order to get the correct anomaly we have to ensure that we are not in the regime of 
the sudden approximation. If we are in the intermediate region we will see a crossover from one limit to 
another limit. So far we have only studied fermions and (abelian) gauge fields in 1+1 dimensions. It was 
pointed out in that the anomaly in 3+1 dimensions factorizes into a 1+1 dimensional part and an extra 
factor coming from the additional dimensions. We briefly review the argument in | p2[ . To get the anomaly 
in 3+1 dimensions we first turn on a magnetic field in the, say z, direction. This leads to the usual Landau 
levels for the fermions which are labelled by integers. We then turn on an electric field paralell to the 
magnetic field. The important point is that the fermions in the lowest Landau level in the presence of this 
electric field behave like fermions in the 1+1 dimensional case that we have just analyzed. Hence the same 
1+1 dimensional anomaly is present but with an additional degeneracy factor coming from the Landau 
levels. As shown in the argument goes through for the lattice Dirac sea case for the case of an uniform 
magnetic field. The degeneracy factor is a geometrical quantity that is in general dependent on the details 
of the lattice Hamiltonian which is more complicated for Ginsparg- Wilson fermions. However, when the 
Ginsparg- Wilson operator is constructed as aDcw = (1 ^ VWa^ with A = 1 — aDw, the degeneracies of 
Dgw and Dw are the same, and in the limit of zero lattice spacing the degeneracy is just 

LiL2gH/{27T) , (47) 

the number of states in the square L1L2 perpendicular to the magnetic field (H). The above factor simply 
multiplies the 1 + 1 dimensional anomaly and gives the correct anomaly in 3+1 dimensions. 

The main aim of this note was to show that a Hamiltonian analysis of Ginsparg- Wilson fermions leads to a 
non-zero rate of chiral charge and gives the anomaly equation in the continuum limit. The doubler modes 
are suppressed by the Wilson parameter, infact the Wilson parameter r plays a crucial role in yielding 
the correct anomaly. A quantum mechanical analysis supplemented by an adiabatic approximation was 
necessary to get the anomaly. It is noteworthy that if we are not in the adiabatic regime we will get 
other contributions (a(t) and higher time derivatives) to the chiral charge and this will not reproduce the 
anomaly equation. It may be useful to compare our derivation with that of the overlap method. In the 
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overlap method the anomaly is extracted by looking at the scalar product of the ground states of two 
different many body Hamiltonians, whereas in our approach we study the dynamical picture behind the 
anomaly by using the properties of the Dirac sea in an external electric field in the adiabatic limit. 



16 



Appendix 

In this appendix we collect together some expressions which are necessary to get the function C{k). 

. /l-2 g (m+ '' '1-^""'° ) +a2 (n^+ - (i-cos(a k ^y ^^2 

m = — ^ — ^ — - — ^ (48) 



9ik) = . -("^-^ (49) 

( 1 — m — ' '^)) j ( — 2 r sin(a k) + 2ar (m + ^ (^~™^("^)) j sin(a k) + 2 cos(a k) sin(a k) ) 

/(fc) = A lA V J__ _ J_ + 

2a(l-2a (m + ^^^i^2££(£M) + (m + + sin(a fc)^^ 

r sm{a k) 



a (^-2a (m+ ''^^-7^°^» ) +a2 (m + + sin(afc)'^ 

— sin(a k) 



2\ 2 



2a(l-2a (m + ^^-7^°^" ) + ^ ra-c£s(afc)) ^ +sin(aA;)' 

cos(a k) 

a{l-2a (m + + a2 (^^ ^ Hi-^o.{ak)) \^ +sin(afc)' 



i;2(fc) 



sin(a /c) 



2 



a2 (^l-2a (m + ^:ii^7^) + a2 (m + ^^ii^22£(£M^% sin(a A;)' 



1 - 



l-2a(m+ ''''-°°°'°'°» )+a2 (nt+ »■ '=» ) %sin(a fc)^ 



C(fc) = (l/(2£;(fc))(5(fc)/(fc) - /(fc)ff(fc)) (50) 
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FIG. 1. C{k) at m=0 and r=0. 



FIG. 2. C{k) at m=0 and r=0.4 



FIG. 3. C{k) at m=0 and r=0.6 



FIG. 4. C{k) at m=0 and r=1.0 



FIG. 5. C{k) at m=0 and r=5.0 



FIG. 6. E{k) at m=0 and r=0.0 



FIG. 7. E{k) at m=0 and r=0.4 



FIG. 8. E{k) at m=0 and r=0.6 



FIG. 9. E{k) at m=0 and r=1.0 



FIG. 10. E{k) at m=0 and r=5.0 



FIG. 11. E{k) at m=5 r=0 



FIG. 12. E{k) at m=5 r=5 
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